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Abstract 

A  nanofluid  is  a  suspension  of  ultrafine  particles  in  a  conventional  base  fluid  which  tremendously 
enhances  the  heat  transfer  characteristics  of  the  original  fluid.  Furthermore,  nanofluids  are  expected 
to  be  ideally  suited  in  practical  applications  as  their  use  incurs  little  or  no  penalty  in  pressure  drop 
because  the  nanoparticles  are  ultrafine,  therefore,  appearing  to  behave  more  like  a  single-phase  fluid 
than  a  solid-liquid  mixture.  About  a  decade  ago,  several  published  articles  focused  on  measuring  and 
determining  the  effective  thermal  conductivity  of  nanofluids,  some  also  evaluated  the  effective 
viscosity.  There  are  only  a  few  published  articles  on  deriving  the  forced  convective  heat  transfer  of 
nanofluids.  The  purpose  of  this  article  is  to  summarize  the  published  subjects  respect  to  the  forced 
convective  heat  transfer  of  the  nanofluids  both  of  experimental  and  numerical  investigation. 
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1.  Introduction 

Improvements  to  make  heat  transfer  equipment  more  energy  efficient  would  need  to 
focus  on  miniaturization  on  the  one  hand  and  an  astronomical  increase  in  heat  flux  on  the 
other.  Heat  transfer  fluids  such  as  water,  mineral  oil  and  ethylene  glycol  play  a  vital  role  in 
many  industrial  processes,  including  power  generation,  chemical  processes,  heating  or 
cooling  processes,  and  microelectronics.  The  poor  heat  transfer  properties  of  these 
common  fluids  compared  to  most  solids  is  a  primary  obstacle  to  the  high  compactness  and 
effectiveness  of  heat  exchangers.  The  essential  initiative  is  to  seek  the  solid  particles  having 
thermal  conductivities  several  hundreds  of  times  higher  than  those  of  conventional  fluids. 
An  innovative  idea  is  to  suspend  ultrafine  solid  particles  in  the  fluid  for  improving  the 
thermal  conductivity  of  a  fluid.  Many  types  of  particle,  such  as  metallic,  non-metallic  and 
polymeric,  can  be  added  into  fluids  to  form  slurries.  However,  the  usual  slurries,  with 
suspended  particles  in  the  order  of  millimeters  or  even  micrometers  may  cause  some  severe 
problems.  The  abrasive  action  of  the  particles  causes  the  clogging  of  flow  channels,  erosion 
of  pipelines  and  their  momentum  transfers  into  an  increase  in  pressure  drop  in  practical 
applications.  Furthermore,  they  often  suffer  from  instability  and  rheological  problems.  In 
particular,  the  particles  tend  to  settle  rapidly.  Thus,  although  the  slurries  give  better 
thermal  conductivities,  they  are  not  practical. 

The  use  of  particles  of  nanometer  dimension  was  first  continuously  studied  by  a  research 
group  at  the  Argonne  National  Laboratory  around  a  decade  ago.  Choi  [1]  was  probably 
the  first  one  who  called  the  fluids  with  particles  of  nanometer  dimensions  ‘nano-fluids’. 
Compared  with  suspended  particles  of  millimeter-or-micrometer  dimensions,  nanofluids 
show  better  stability  and  rheological  properties,  dramatically  higher  thermal  conductiv¬ 
ities,  and  no  penalty  in  pressure  drop.  Several  published  literature  have  mainly  focused  on 
the  prediction  and  measurement  techniques  in  order  to  evaluate  the  thermal  conductivity 
of  nanofluid.  Recently,  Trisaksri  and  Wongwises  [2]  reviewed  the  literature  on  the  general 
heat  transfer  characteristics  of  nanofluids.  It  was  found  that  only  a  few  papers  have 
discussed  on  the  convective  heat  transfer  of  nanofluids.  The  purpose  of  this  article  is  to 
review  the  literature  mentioning  convective  heat  transfer  of  nanofluids,  including  the 
experimental  and  theoretical  investigations. 

2.  Experimental  approach 

A  several  years  ago,  the  nanofluid  has  been  found  to  be  an  attractive  heat  transport 
fluids.  It  has  exhibited  a  significant  potential  for  heat  transfer  augmentation  relative  to  the 
conventional  pure  fluids.  It  has  been  expected  to  be  suitable  for  the  engineering  application 
without  severe  problems  in  pipeline  and  with  little  or  no  penalty  in  pressure  drop. 
However,  the  earlier  researches  revealed  that  most  studies  of  the  transport  properties  of  the 
nanofluid  were  focused  on  the  thermal  conductivity.  The  convective  heat  transfer  of 
nanofluids  has  received  comparatively  little  attention  in  literature.  There  are  few 
publications  dealing  with  the  convective  heat  transfer  of  nanofluids.  The  most  productive 
studies  have  been  continuously  carried  out  by  the  following  researchers. 
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Nomenclature 

A  cross  section  area,  m2 

C  concentration  of  liquid  solution,  kgw/kgsoi 

Cp  specific  heat,  J/kg  K 

d  nanoparitcle  diameter,  m 

D  inner  diameter  of  the  tube,  m 

g  acceleration  of  gravity,  m/s2 

h  heat  transfer  coefficient,  W/m2  K 

H  humidity  ratio  of  the  air,  kgw/kga 

k  thermal  conductivity  of  the  nanofluid,  W/m  K 

L  length  of  the  test  tube,  m 

m  mass  flow  rate,  kg/s 

Nu  Nusselt  number 

A P  pressure  drop,  Pa 

Pe  Peclet  number 

Pr  Prandtl  number 

q  heat  flux,  W/m2 

Q  heat  transfer  rate,  W 

r  radius  of  tube,  m 

R  radius  of  disc,  m 

Re  Reynolds  number 

S  perimeter,  m 

/  thickness,  m  or  time  in  the  Lattice  Boltzmann  method 

T  temperature,  °C 

u  mean  velocity,  m/s 

U  overall  heat  transfer  coefficient,  W/m2  K 

v  velocity  in  y-direction,  m/s 

V  volume,  m3 

w  velocity  in  z-direction,  m/s 

x  distance  along  the  tube,  m  or  x-coordinate 

y  v-coordinate 

z  z-coordinate 

a  thermal  diffusivity,  m2/s 

ft  thermal  dispersion  coefficient,  N/m2  K 

£  rate  of  dissipation  (for  k—e  turbulent  model) 

X  friction  factor 

co  tangential  coordinate,  rad 

cj)  volume  fraction 

< P  inclination  angle,  degree 

t/y  thermal  diffusivity  coefficient 

W  mass  diffusivity,  m2/s 

k  turbulent  kinetic  energy  (for  k-s  turbulent  model) 
p  density,  kg/m3 

t  tensor  of  the  dispersed  thermal  conductivity 

v  viscosity,  m2/s 
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Subscript 

a  annulus  or  air 

d  particle 

des  desiccant 

D  dispersion  effect 
f  fluid 

hf  heating  fluid 

i  inside 

in  inlet 

lrn  log  mean 

m  mean 

o  out 

out  outlet 

nf  nanofluid 

r  radial  direction 

T  turbulent 

w  tube  wall 

x  axial  direction 


Li  and  Xuan  [3],  and  Xuan  and  Li  [4]  presented  an  experimental  system  to  investigate 
the  convective  heat  transfer  coefficient  and  friction  factor  of  nanofluids  for  laminar  and 
turbulent  flow  in  a  tube.  The  sample  nanofluid  used  in  this  study  was  deionized  water  with 
a  dispersion  of  Cu  particles  with  below  lOOnm  diameter.  The  nanofluids,  with  different 
particle  volume  percentages,  were  used  to  determine  the  effect  of  the  nanoparticle 
concentration  on  the  heat  transfer  coefficient,  the  percentages  being  0.3,  0.5,  0.8,  1,  1.5  and 
2%.  The  Reynolds  number  of  the  nanofluids  varied  within  the  range  of  800-25,000.  The 
experimental  apparatus  was  mainly  composed  of  a  reservoir  tank,  a  pump,  a  pipeline,  a 
test  section,  a  cooler,  and  a  fluid  collection  tank.  Test  section  was  a  straight  brass  tube  with 
10  mm  inner  diameter  and  800  mm  length  and  a  3.5  kW  electric  heater  was  used  to  obtain  a 
constant  wall  heat  flux  boundary  condition.  To  prevent  the  aggregation  of  the 
nanoparticles  in  the  base  fluid,  a  fatty  acid  was  used  as  a  dispersant  to  cover  the 
nanoparticles.  The  experimental  apparatus  of  their  work  was  depicted  in  Fig.  1. 

The  heat  transfer  coefficient  of  the  nanofluids  is  calculated  from  the  following  equations 

Nunf  =  h'f  (1) 

E-  nf 


^nf  — 


g 

Fw  -  Tf 


(2) 


The  experimental  results  illustrated  that  the  convective  heat  transfer  coefficient  of  the 
nanofluids  varied  with  the  flow  velocity  and  volume  fraction  as  well  as  being  higher  than 
the  base  fluid  at  the  same  conditions.  Compared  with  water,  the  Nusselt  number  of  the 
nanofluid  with  2%  volume  fraction  of  Cu  particles  was  60%  higher.  Furthermore, 
the  Dittus-Boelter  equation  was  not  valid  for  the  prediction  of  the  Nusselt  number  of  the 
nanofluids  at  various  volume  fractions.  From  the  experimental  data,  new  heat  transfer 


W.  Daung thong suk,  S.  Wongwises  /  Renewable  and  Sustainable  Energy  Reviews  11  (2007)  797-817  801 


DC  Power  Supply 


Fig.  1.  The  experimental  apparatus  of  Xuan  and  Li  [3]  and,  Li  and  Xuan  [4]. 


correlations  for  the  prediction  of  the  heat  transfer  coefficient  of  nanofluids  flowing  in  a 
tube  could  be  established  in  the  following  form: 

For  laminar  flow 


Nun  f  =  0.4328(1.0+  11. 285/'754F4218)Re“f333Fr”f4 


For  turbulent  flow 


Nun  f  =  0.0059(1.0  +  7.6286^)' 


0.6886^0.001 


)Re' 


0.9238p;.0.4 


(3) 


(4) 


where  Nu  is  the  Nusselt  number,  Pe  is  the  Peclet  number.  Re  is  the  Reynolds  number,  and 
Pr  is  the  Prandtl  number. 


PeA  = 

0Cnf 

umD 

R^nf  — 


Pnf 


and 


Prni  =  — 

CXnf 


(5) 


(6) 


(7) 


In  order  to  calculate  these  parameters,  the  thermal  diffusivity  of  the  nanofluid  was 
defined  as 


0Cnf 


^nf 


^nf 


{pCp)  nf  O-^)(pCp)f  +  0(pCp)d 


(8) 


In  the  case  of  the  pressure  drop,  the  results  indicate  that  the  friction  factor  of  the 
nanofluid  was  almost  equal  to  those  of  water  under  the  same  flow  velocity  and  did  not  vary 
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with  the  volume  fraction.  This  implies  that  the  nanofluids  cause  no  penalty  to  pump  power 
and  the  single-phase  friction  factor  correlation  can  be  used  for  more  accurate  predictions. 
The  friction  factor  of  the  nanofluids  was  determined  from  the  following  equation: 


Anf  = 


A Pn(D  2g 

L  < 


(9) 


Yang  et  al.  [5]  reported  experimental  results  which  illustrated  the  convective  heat 
transfer  coefficient  of  graphite  nanoparticles  dispersed  in  liquid  for  laminar  flow  in  a 
horizontal  tube  heat  exchanger.  Nanoparticles  of  different  sources  were  used  in  this  study 
and  focused  on  an  aspect  ratio  (1/d)  of  about  0.02  (like  a  disc),  because  the  addition  of  large 
aspect  ratio  particles  into  a  fluid  may  dramatically  increase  viscosity  compared  with  the 
single-phase.  In  this  study,  two  series  of  nanofluids  with  different  base  fluids  were  used. 
The  experiments  were  performed  within  the  following  ranges:  the  volume  flow  rates  were 
between  62-507  cm3/min,  the  Reynolds  number  varied  between  5-110,  and  the  heating 
fluid  temperature  between  50-70  °C.  This  article  intended  to  study  the  effects  of  the 
Reynolds  number,  volume  fraction,  temperature,  nanoparticle  source,  and  type  of  base 
fluid  on  the  convective  heat  transfer  coefficient.  The  experimental  system  used  in  this  study 
was  shown  schematically  in  Fig.  2. 

The  following  calculation  procedure  was  used  to  evaluate  the  heat  transfer  coefficient  of 
the  nanofluids. 

The  heat  transfer  rate  into  the  nanofluids  is 


(2nf  =  mCp(Tout  -  T in)  (10) 

The  heat  transfer  from  the  heating  fluids  is 

Qhf  —  mhfCphf(T0Ut  —  7’in)hr  (11) 

The  energy  balance  between  the  heating  fluid  and  the  process  fluid  were  differences 
around  25%.  This  is  common  for  a  small-scale  testing  equipment. 


T  T 

A  in  ±  out 


Fig.  2.  Schematic  diagram  of  the  experimental  apparatus  of  Yang  et  al.  [5], 
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The  overall  heat  transfer  resistance  relationship  is  used  to  determine  the  heat  transfer 
coefficient  of  the  nanofluids,  hn{: 

1  1  |  A.,  A,  |  1 

U~  hn[(Ai/A0)  +  2kmDl+h0 

The  overall  heat  transfer  coefficient  can  then  be  determined  from 


£/  = 


Q 

A0ST\m 


(13) 


The  Monrad  and  Pelton’s  equation  for  turbulent  flow  in  annuli  was  used  to  calculate  the 
outside  heat  transfer  coefficient,  /;G,  as  follows: 


Nu  =  0.020Rea8Pr1/3  (Z>2/Z>i)a53  ( 14) 

Subscript  2  and  1  mean  the  outer  and  inner  diameter  of  the  annulus,  respectively.  The 
equivalent  diameter,  De,  was  used  instead  of  the  typical  diameter,  D. 


De  = 


Dr 


Dr 


D0 


(15) 


The  experimental  results  illustrated  that  the  heat  transfer  coefficient  increased  with  the 
Reynolds  number  and  the  particle  volume  fraction.  For  example,  at  2wt%  the  heat 
transfer  coefficient  of  the  nanofluids  was  moderately  increased  compared  with  the  base 
fluid,  but  at  2.5  wt%  it  was  22%  higher  than  the  base  fluid  at  50  °C  fluid  temperature  and 
15%  higher  at  70  °C.  Several  factors  may  lead  to  smaller  improvements  in  heat  transfer 
enhancement  of  nanofluids  with  increased  fluid  temperature.  One  is  the  rapid  alignment  of 
nanoparticles  in  less  viscous  fluids,  leading  to  less  contact  between  the  nanoparticles.  A 
second  factor  is  that  depletion  of  particles  in  the  near-wall  fluid  phase  could  lead  to  a  lower 
thermal  conductivity  layer  near  the  wall.  Two  graphite  nanoparticle  sources  at  same 
particle  loading  give  different  heat  transfer  coefficient,  as  in  the  case  of  different  base 
fluids.  This  may  have  been  caused  by  particle  shape,  morphology  or  surface  treatment. 
This  means  that  one  type  of  nanoparticle  is  more  effective  than  another  in  augmenting  the 
heat  transfer  coefficient  of  the  nanofluid.  Moreover,  the  correlation  established  by  Li  and 
Xuan  [4]  for  laminar  flow  of  nanofluids  gave  a  higher  heat  transfer  coefficient  than  that 
calculated  from  the  experimental  data.  This  may  be  because  the  effect  of  particle  aspect 
ratio  was  not  included.  Similarly,  the  single-phase  correlation  of  laminar  flow  in  tubes 
failed  to  predict  the  convective  heat  transfer  coefficient  of  nanofluids. 

All  of  the  experimental  data  were  used  to  develop  a  new  heat  transfer  correlation  for  the 
prediction  the  heat  transfer  coefficient  of  laminar  flow  nanofluids  in  a  more  convenient 
form  by  modifying  the  Seider-Tate  equation  [6]  as  follow: 


Q  =  Nun[Pr~fl/3 


=  aRc 


b 

nf 


(16) 


where  a  and  b  are  the  various  constants  which  depend  on  the  experimental  conditions. 

Wen  and  Ding  [7]  presented  an  experimental  study  which  evaluated  the  convective  heat 
transfer  coefficient  of  7AI2O3  nanoparticles  suspended  in  deionized  water  for  laminar  flow 
in  a  copper  tube,  focusing  in  particular  on  the  entrance  region.  The  various  concentrations 
of  nanoparticles  (0.6,  1.0,  and  1.6%)  were  tested  under  a  constant  wall  heat  flux.  The 
experimental  apparatus  mainly  consisted  of  a  straight  copper  tube  with  970  mm  length. 
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4.5  +  0.02  mm  inner  diameter,  and  6.4  +  0.05  mm  outer  diameter  as  a  test  section.  A  300  W 
silicon  rubber  flexible  heater  was  used  to  keep  a  constant  heat  flux  boundary  condition.  A 
peristaltic  pump  was  used  to  deliver  a  maximum  flow  rate  of  101/min.  The  flow  rate  was 
controlled  by  adjusting  rotational  speed  of  the  pump.  In  order  to  stabilize  the 
nanoparticles,  sodium  dodecylbenzene  sulfonate  (SDBS)  was  used  as  a  dispersant.  Their 
experimental  apparatus  was  shown  schematically  in  Fig.  3. 

The  heat  transfer  coefficient  and  the  Nusselt  number  of  the  nanofluids  were  calculated 
from  the  following  equations: 


N«nf,x  = 


kn( 


(17) 


/inf,x  =  „  q  Tf  (18) 

r  w,x  —  1  f,x 

The  energy  balance  was  used  to  determine  the  fluid  temperature  profile  along  the  test 
tube  as  follows: 


T  fjX  —  T  in  + 


qSx 

pCpuA 


(19) 


The  results  showed  that  the  local  heat  transfer  coefficient  varied  with  the  Reynolds 
number  and  volume  fraction.  In  particular,  the  entrance  region  illustrated  that  the  use  of 
nanofluids  affected  a  pronounced  increase  in  the  heat  transfer  coefficient  causing  a 
decrease  in  the  thermal  boundary  layer  thickness  and  will  decrease  along  the  tube  length. 


DC  Power  Supply 


Fig.  3.  Schematic  diagram  of  the  experimental  apparatus  of  Wen  and  Ding  [7]. 
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This  behavior  implied  that  it  might  be  possible  to  create  a  ‘smart  entrance’  region  to  meet 
the  highest  performance  of  nanofluids.  Furthermore,  the  Shah  correlation  for  laminar  flow 
and  the  Dittus-Boelter  equation  for  turbulent  flow  do  not  coincide  with  the  Nusselt 
number  between  the  experimental  results  and  the  calculated  values. 


3.  Mathematical  approach 


Similar  to  the  experimental  investigation,  many  researchers  devoted  to  formulate  the 
appropriate  correlations  for  prediction  of  the  thermal  conductivity  of  the  nanofluid.  It  is 
difficult  to  establish  any  formulated  theory  for  the  prediction  of  the  accurate  heat  transfer 
behavior  of  the  nanofluid  because  the  several  complex  phenomena  may  coexist  in  the  main 
flow  e.g.  sedimentation,  Brownian  force,  dispersion,  and  friction  force  between  the  fluid 
and  the  solid  particles.  Some  researchers  treat  the  nanofluid  more  like  a  fluid  rather  than  a 
liquid-solid  mixture.  There  are  few  articles  taking  the  nanofluid  into  account  as  the  multi¬ 
phase  feature.  Some  detailed  reviews  mentioned  above  are  shown  as  follows: 

Xuan  and  Li  [8]  presented  a  method  for  preparing  some  nanofluid  samples  and 
introduced  a  theoretical  study  of  thermal  conductivity.  They  used  the  transient  hot-wire 
method  for  measurements.  They  then  suggested  a  dispersion  model  or  single-phase  model 
to  describe  the  heat  transfer  performance  of  nanofluids  flowing  in  a  tube.  The 
nanoparticles  used  in  this  article  were  Cu  particles  with  lOOnm  diameter,  respectively, 
suspended  in  water  and  mineral  oil.  Oleic  acid  and  laurate  salt  were  used  as  dispersants  in 
order  to  stabilize  the  particles  in  the  mineral  oil  and  water,  respectively.  For  Cu-mineral 
oil,  the  results  showed  that  the  suspension  reached  a  stabilized  state  within  about  1  week 
and  had  no  sedimentation  with  a  fill  of  22  wt%.  In  the  case  of  Cu-water  nanofluid 
suspension  with  9  wt%  laurate  salt,  vibrated  by  a  ultrasonic  vibrator,  the  results  indicated 
that  the  nanoparticles  can  suspend  more  30  h  than  in  stationary  state.  This  behavior 
implied  that  the  fluid  viscosity  might  be  an  important  factor  affecting  the  dispersion  and 
stability  of  the  particles.  Moreover,  this  study  illustrated  that  the  Hamilton  and  Crosser 
model  [9]  and  Wasp  model  [10]  for  predicting  thermal  conductivity  can  be  applied  to 
obtain  an  approximate  solution  in  the  absence  of  a  sophisticated  formula.  Using  the 
transient  hotwire  method  to  measure  the  thermal  conductivity,  the  results  showed  that  the 
thermal  conductivity  remarkably  varied  with  the  volume  fraction  of  the  particles.  There 
are  two  different  approaches  to  analyze  the  heat  transfer  enhancement.  One  is  the  two- 
phase  model  and  other  is  the  single-phase  model.  Although  the  first  model  describes  an 
understanding  of  the  functions  of  both  liquid  phase  and  solid  phase  in  the  heat  transfer 
process,  it  needs  a  long  time  for  computation  and  a  high  performance  computer.  The 
second  model,  assuming  that  both  the  liquid  and  particle  phases  are  in  thermal  equilibrium 
and  flow  at  the  same  velocity,  is  simpler  and  needs  less  computation  time.  The  purpose  of 
their  work  aimed  to  develop  a  modified  single-phase  model  to  describe  the  heat  transfer 
process  of  nanofluids  flowing  in  tubes.  The  governing  equations  were  expressed  as  follow: 


dT  dT 

— - h  U  — —  — 

dt  dx 


^_\^T  1 1 


(pCp)nf/  5x2  r8r 


<t>  + 


Pr 


dt 


(pCpW  dr 


(20) 


If  the  axial  temperature  distribution  is  ignored,  Eq.  (20)  can  be  reduced  as 


dT  dT  _  Id 
dt  U  dx  r  dr 


(pCp)J 


r 


dT' 

~d7 


(21) 
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The  boundary  conditions  for  this  situation  are  as  follows: 

T\r=R  =  Tw  and  T\x=q  =  T0  (22) 

The  application  of  the  separation  of  variables  method  for  plug  flow,  the  following 
solution  was  obtained: 


T  -  77 


n=  1 


^-O^x/Pe  Jo(®mr) 
Me  m)$m 


(23) 


where  r  —  r/R ,  x  —  x/L,  Nu  =  2Rh//c*f,  Pe*  =  uL/<ft*,  Pe  =  Pe*(R/L,y  and  eigenvalues 
0m’s  are  the  positive  roots  of  the  following  equation: 


J  0  (dm)  —  0 

The  heat  balance  can  be  expresses  as 


r=R 


=  h(Th  -  rw) 


The  analytical  correlation  of  the  Nusselt  number  is  then  given  by 


(24) 

(25) 


Nu  = 


Em=le 


-8}„x/Pe 


Eoo 

,  e  n 

m=  1 


vn/qi 


(26) 


For  presence  of  the  dispersion  in  axial  and  radial  direction,  Eq.  (20)  and  the  modified 
Dankwerts’s  boundary  conditions  (Eqs.  (27a)  and  (27b))  in  the  axial  direction  and  the 
boundary  condition  in  radial  direction  Eq.  (22)  give  solution  (Eq.  (28)). 


—k 


* 

nf 


dT 

dx 


—  uApCp(TQ  —  T)  at  x  —  0, 


(27a) 


— —  =  0  at  x  —  L 

dx 


(27b) 


T  -  7\v  _  2  f-  x(x)  Me„f) 
To-Tw  -  2^  [X(0)  -  r (0 )/Pe*}  J x(6m)em 

where  X(x)  and  X'(x)  are  defined  as  follows: 

A  (  v)  =  rn2em2+mi*  -  w1emi+m2je 


X'{x)  =  m\rn2(em2+mi*  -  emi+m2*) 

Pe*  ±  J  Pe*2  +  A92m(L/R)2 
»h,2  = - j - 

The  eigenvalues  are  described  in  Eq.  (24).  Thus,  the  Nusselt  number  correlation  for  the 
nanofluids  is  obtained  from 

Nl[  _  [ym  -  x\o)/Pe*] 

Em=i  [n*)l  [3f(0)  -  X'(0 )/Pe*\ }  /e2m 


(29) 
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It  should  be  noted  that  Eqs.  (26)  and  (29)  are  derived  for  constant  wall  temperature  and 
need  the  iteration  procedure  to  solve  the  problem.  The  Nusselt  number  for  the  other 
boundary  conditions  can  derive  with  the  same  method. 

Finally,  this  study  expressed  that  either  theoretical  expressions  or  experimental  work 
needed  to  explain  the  heat  transfer  process  of  nanofluids,  in  which  several  phenomena  may 
coexist  in  the  main  flow. 

Xuan  and  Roetzel  [11]  presented  two  different  approaches  to  determine  some 
fundamentals  to  predict  the  convective  heat  transfer  coefficient  of  nanofluids  under 
the  assumption  that  they  behaved  like  a  single-phase  liquid  rather  than  a  normal 
solid-liquid  mixture.  They  included  both  the  effects  of  transport  properties  and  thermal 
dispersion  of  the  nanofluid.  The  first  model  treated  the  nanofluid  as  a  single-phase 
fluid  and  the  second  model  as  a  multi-phase  and  dispersed  fluid.  In  a  conventional 
approach,  the  assumptions  were  no  slip  velocity  between  the  particles  and  liquid  phase 
and  that  they  were  in  a  thermal  equilibrium  state.  The  nanofluid  behaved  like  a  common 
pure  liquid.  This  implied  that  the  whole  equation  of  continuity,  motion  and  energy  for  a 
single-phase  fluid  had  been  applied  directly  to  the  nanofluid.  Similar  to  the  Xuan  and  Li 
model  [8],  the  Hamilton  and  Crosser  model  [9]  and  Wasp  model  [10]  for  predicting  the 
thermal  conductivity  of  nanofluid  give  rough  estimations.  For  a  modified  conventional 
approach,  called  the  dispersion  model,  the  slip  velocity  between  the  solid-liquid  mixture 
may  not  be  zero  causing  many  factors,  such  as  Brownian  force,  gravity,  friction  force 
between  the  solid-liquid  mixture,  Brownian  diffusion,  sedimentation,  and  dispersion  to 
happen  simultaneously  in  the  main  flow.  The  details  of  two  approaches  are  expressed 
as  follows: 

For  the  conventional  approach,  the  general  form  of  the  energy  equation  for  the 
incompressible  pure  fluid  with  negligible  viscous  dissipation  is 


(30) 


This  equation  is  suitable  to  explain  the  heat  transfer  process  of  the  nanofluid.  This 
means  that  the  heat  transfer  equation  for  the  single-phase  fluids  can  extend  to  the 
nanofluid  directly.  It  should  be  noted  that  the  thermal  properties  used  in  this  equation  are 
those  of  nanofluids.  For  this  reason,  the  dimensionless  heat  transfer  correlations  of  the 
single-phase  fluid  can  be  applicable  for  the  nanofluid,  so 

Nu  =  3.66  (For  fully  developed  laminar  flow  in  tube,  Tw  =  const.)  (31a) 

Nu  =  0.023Re°'8Pr1/3  (For  turbulent  flow  in  tube)  (31b) 

For  the  modified  conventional  approach,  the  irregular  movement  of  the  ultrafine 
particles  led  to  little  disturbances  of  both  the  temperature,  T ,  and  velocity  of  the 
nanofluid,  u' .  Hence,  the  real  temperature  and  velocity  averages  are  expressed  as 


T  =  {Tf  +  T 


(32a) 


u  =  („}/  +  u' 


where  (T/  =  (1/Ff)  [Vf  TdV,  (u)f  =  (1/Ff)  jV(  udV  and  (1/Ff)  jVf  T'dV  =  0 


(32b) 
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By  the  procedure  demonstrated  by  Kaviany  [12]  and  assuming  that  the  boundary  surface 
between  the  fluid  and  the  particles  is  so  small,  Eq.  (30)  can  be  rewritten  as 


(pC  p)nf 


'8{T)f 

dt 


+  (h/V  ■  (Tf 


V-  (kn(V{Tf)  -  (pCp)nfV(ii,r>/ 


(33) 


The  chaotic  movement  of  the  ultrafine  particles  in  the  flow  gives  the  thermal  dispersion 
(the  second  term  on  the  right  hand  side  of  Eq.  (33)).  Therefore,  the  heat  flux  resulted  from 
the  dispersion  in  the  nanofluid  flow  can  be  expressed  as: 

qn  =  (pCP)nf(«/r)/  =  -ID  .  V(2y  (34) 

where  t  is  the  tensor  of  the  dispersed  thermal  conductivity.  Thus,  Eq.  (33)  can  be 
rewritten  as 


d{T)' 


dt 


+  (u/V-(T)f 


=  V- 


Then,  the  energy  balance  equation  for  the  nanofluid  can  be  given  as 

8{T)' 


(pc  p)nr 


8t 


+  (pC P)nf(«>  •  V(T)  =  (/0CP)nfV(.A  ■  V(T» 


(35) 


(36) 


where  i p  may  be  called  the  total  effective  thermal  diffusivity  tensor.  It  combines  the 
molecular  effect  and  the  thermal  dispersion  effect,  i//  is  defined  as  follows: 


<A  -  ,  +  (1  -<A)<Ad 

(pCp)nf 


(37) 


It  is  anticipated  that  the  thermal  dispersion  tensor  is  a  function  of  the  flow  pattern, 
properties  of  the  base  fluid  and  the  nanoparticles,  dimensions  and  shape  of  the 
nanoparticles,  and  volume  fraction  of  the  nanoparticles. 

The  energy  equation  for  the  nanofluid  flowing  inside  a  tube  for  this  situation  may  be 
expressed  as 


8T  8T  _  1  8 
dt  +  U  8x  r  dr 


0Cnf  ' 


cD,r 


-VS 


(pCp)n{J  dr  _ 

For  taking  the  thermal  dispersion  in  the  axial  direction  into  account,  so 


8T  8T  _  1  8 
dt  dx  r  dr 


0inf  ' 


fep.r  \r8T 
(pCp)nf/  dr  _ 


+ 


dx 


0Cnf  ■ 


k 


D,jc 


8T 


(, pCv)ntJ  dx 


(38) 


(39) 


The  thermal  diffusivity  coefficients  can  be  found  from  the  derivation  of  Taylor  [13]  and 
Aris  [14].  Unfortunately,  there  is  no  theoretical  and  experimental  work  established  on  the 
dispersed  thermal  diffusivity  of  the  nanofluid.  The  kD  appearing  in  Eqs.  (38)  and  (39)  can 
be  calculated  by  Plumb  [15]  and  Elunt  and  Tien  [16]. 

To  take  the  effect  of  thermal  dispersion  into  account,  the  surface  heat  flux  can  be 
expressed  as 


8T 

dr 


I  r=R 


(40) 


Similar  to  the  procedure  proposed  by  Nield  and  Bejan  [17],  the  dimensional  analysis  of 
the  above  equation  formally  resulted  in  the  heat  transfer  correlation  of  the  nanofluid  as 
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follows: 


Nux  =  [1  +  C*Pef'(0)]9'(0)Rem 


(41a) 


with  considering  the  flow  pattern,  Eq.  (41a)  can  be  modified  to  Eq.  (41b) 
Nux  =  [1  +  C*  Pe'f  (0)]  0'(  0 )  Re'" 


(41b) 


In  which  Pe  =  (rum/ Onf),  C*  is  an  unknown  constant  determined  from  experiment,  m  and 
n  depend  on  the  flow  pattern,  /  is  dimensionless  velocity,  0  is  dimensionless  temperature, 
and  f  and  O'  are  the  derivatives  of  the  dimensionless  velocity  and  dimensionless 
temperature,  respectively. 

Maiga  et  al.  [18]  presented  a  mathematical  formulation  and  numerical  method  to 
determine  the  forced  convective  heat  transfer  and  wall  shear  stress  for  the  laminar  and 
turbulent  regions  of  water-yAl203  and  ethylene  glycol-}’ AFO3  flowing  inside  a  uniformly 
heated  tube.  The  solid-liquid  mixture  takes  a  single-phase  behavior  into  account,  so  the 
slip  velocity  between  the  phases  was  neglected.  Furthermore,  the  local  thermal  equilibrium 
of  the  mixture  and  symmetry  in  flow  were  considered.  For  the  turbulent  flow  region,  the 
Reynolds-averaged  Navier-Stokes  equation  and  K-e  turbulent  model  were  used  to  describe 
the  stresses  and  heat  flux  of  the  nanofluids.  The  conditions  of  the  numerical  simulations 
were  as  follow:  the  heated  tube  had  0.01  m  diameter  and  1.0  m  length.  In  the  situation  of 
laminar  flow,  the  Reynolds  number  was  fixed  at  250  and  the  constant  heat  flux  was  varied 
between  10-250  W/m2.  For  turbulent  flow,  the  constant  heat  flux  was  fixed  at  500,000  W/ 
m2  and  the  Reynolds  number  was  varied  in  the  range  of  10,000-50,000.  The  problem 
configuration  for  this  study  is  shown  in  Fig.  4. 

Under  the  above  conditions  the  governing  equations  depicted  in  a  ( R ,  co  and  Z) 
coordinate  system  are  as  follows: 


V  ■  {pit)  =  0 

V  ■  (puiij)  =  —dP/dXj  +  V  ■  (/(Vm,)  +  St 

V  ■  (puCpT)  =  V  ■  (/c V T) 


(42) 

(43) 

(44) 


where  u  =  (uR,  u,,„  uz)  is  the  velocity  vector,  Xt  represents  a  spatial  direction  (X,-  =  R,  co 
and  Z),  and  St  represents  the  remaining  viscous  term  corresponding  to  the  coordinate 
system. 

—  for  the  radial  direction,  R(i=  1); 


(45) 


x 


Fig.  4.  The  problem  configuration  in  the  study  of  Maiga  et  al.  [18]. 
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for  the  tangential  direction,  co  ( i  =  2); 
2  ditR 

R 2  dw 


Si  =  p 


7F 


R 


(46) 


—  for  the  axial  direction,  Z  (i  =  3); 

S3  =  0  (47) 

For  the  turbulent  flow  regime,  this  article  shows  the  application  of  the  famous  semi- 
empirical  k— e  turbulent  model.  The  equations  for  this  case  are  described  as  follow: 


V  ■  ( puK )  =  V  ■ 


Wt 

+  K 


+  Gk  -  pe 


(48) 


V  ■  {pm)  —  V  ■ 


P>T,e 


+  Cu{e/k)Gk  +  C2epe2/k 


(49) 


where  GK  is  the  generation  of  turbulence  kinetic  energy  due  to  the  mean  velocity  gradients, 
and  pT  is  the  turbulent  viscosity  calculated  from 


j uT  =  pCliK1/e  (50) 

Various  constants  in  the  above  equations  are  as  follows: 

Cu  =  1-44;  C2e  =  1.92;  C„  =  0.09;  PrT>K  =  1.0;  and  PrT,s  =  1.3  (51) 


The  numerical  results  illustrated  that  the  heat  transfer  and  the  wall  friction  of  nanofluids 
increased  with  increasing  the  particle  concentration  and  the  ethylene  glycol-yAl203  gave 
higher  heat  transfer  enhancement  than  the  water-yAl203.  Then,  the  heat  transfer 
performance  of  the  nanofluids  was  more  pronounced  with  the  increase  of  the  Reynolds 
number  for  the  turbulent  flow  region. 

Roy  et  al.  [19]  proposed  a  numerical  investigation  to  evaluate  the  heat  transfer  and  wall 
shear  stress  for  radial  laminar  flow  in  a  cooling  system  of  water-yAl203  nanofluids 
compared  with  water,  glycol,  and  oil  as  base  fluids.  All  assumptions  used  in  this  study  were 
similar  to  Maiga  et  al.  [18]  who  took  the  nanofluid  as  an  incompressible  fluid.  The 
geometric  configuration  of  this  work  is  shown  in  Fig.  5. 


Outlet 


Z 


Fig.  5.  The  geometric  configuration  in  the  study  of  Roy  et  al.  [19]. 


W  Daung thong suk,  S.  Wongwises  /  Renewable  and  Sustainable  Energy  Reviews  11  (2007)  797-817  811 


The  governing  equations  in  the  cylindrical  coordinate  (R,  Z )  and  symmetry  systems  can 
be  given  as  follows: 

V  ■  (pu)  =  0  (52) 

V  ■  ( puuj )  =  —VP  +  V  ■  (pViij)  +  Si,  i  =  1,2  (53) 

V  -(puCpT)  =  V-(kVT)  (54) 

—  for  the  radial  direction,  R(i  =  1 ); 


Si 


=  ~P 


Ur 

R 2 


—  for  the  axial  direction,  Z(i  =  1); 


(55) 


S2  =  0  (56) 

The  results  show  that  the  heat  transfer  rate  and  wall  shear  stress  were  higher  with 
increasing  particle  concentration  and  Reynolds  number. 

Nguyen  et  al.  [20]  introduced  a  numerical  simulation  to  determine  the  efficiency  of 
water  yAFCF  and  ethylene  glycol  -yAhO,  nanofluids  for  the  cooling  of  a  high-heat  output 
microprocessor  under  laminar  forced  flow  inside  a  heat  sink.  The  simulation  set-up 
consisted  of  a  50  x  50  x  10  mm  rectangular  slot  with  a  3  x  48  mm  fluid  flow  cross-section. 
Assumptions  used  in  this  article  were  laminar  flow,  uniform  velocity  and  temperature 
profile  at  the  inlet  section.  The  contact  area  for  heat  exchange  was  10  x  10  mm.  The 
problem  configuration  for  this  study  was  sketched  in  Fig.  6. 

The  numerical  results  indicate  that  the  average  heat  transfer  coefficient  of  nanofluids 
varied  higher  with  the  volume  fraction  and  Reynolds  number  and  the  ethylene 
glycol-yAl203  nanofluid  gave  higher  heat  transfer  coefficient  than  the  water-yAl203 
nanofluid. 

Ali  et  al.  [21]  reported  the  mathematical  and  numerical  formulation  to  determine  heat 
and  mass  transfer  between  air  and  falling  solution  film  in  a  cross-flow  heat  exchanger  for 
the  dehumidification  and  cooling  process.  The  Cu  nanoparticles  were  added  into  the 
solution  in  order  to  evaluate  their  effect  in  the  augmentation  of  the  heat  and  mass  transfer 
process.  Two  models,  proposed  by  Xuan  and  Roetzel  [11],  were  used  in  this  study. 
Furthermore,  the  assumptions  of  this  article  were  as  follow:  (1)  laminar  and  steady  state 
flow;  (2)  thermal  properties  of  the  air  and  the  solution  are  constant  except  for  the  thermal 


Contact  Area  with  .  Outlet 

Microprocessor  (10x10  mm)  yyr 


Cross-section 


Fig.  6.  The  problem  configuration  in  the  study  of  Nguyen  et  al.  [20]. 
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Flow  of  Liquid 
Film  Desiccant 


conductivity  of  the  solution;  (3)  gravitational  force  of  the  air  is  neglected;  (4)  constant  film 
thickness  and  (5)  fully  developed  velocity  profile  for  fluids.  The  schematic  diagram  of  this 
work  was  shown  in  Fig.  7. 

The  governing  conservation  equations  of  mass,  momentum,  energy,  and  mass  diffusion 
are  illustrated  in  two  portions;  air  and  solution  (liquid  film)  as  follows: 

—  For  the  air 


dwz 

~dZ 

dP a 
dZ 


=  0 


Ha 


p.A C Dn wa  — —  =  ka—^~ 

Pa  Pa  dz  dy2 


,SHa_ 

Wa  dz  dyl 


—  For  the  solution  (liquid  film) 

C^dcs 


dx 


=  0 


Pdesd  4“  hdcs 


d  Wdes 

w 


des 


=  0 


(pCpw)des 


dT  des 
dx 


d 

^des 


(^nf  ^dis) 


dudes' 
^des . 


Wdes 


SCdes 

dx 


f'des 


S^Cdes 

d}’L 


(57) 

(58) 

(59) 

(60) 

(61) 

(62) 

(63) 


(64) 
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The  energy  balance  equation  at  interface  between  air  and  film  desiccant  is 


"b  Pa^a^fg 


djh 


—  ^de: 


STdes 

^.Vdes 


at  ya  =  k,  ydes  =  ?des,  0<x<H  and  0 <z<L 
The  mass  balance  at  the  interface  is 


(65) 


Paf'a 


d_H_ a 

3ja 


=  ~ipV)  des 


5Cdes 

SVdes 


(66) 


at  ya  =  k,  Vdes  =  ?des,  0<x<H  and  ()<:<L 

With  the  suitable  boundary  conditions,  the  velocities  for  both  air  and  film  desiccant  are 
obtained  analytically  from  Eqs.  (58)  and  (62)  as  follows: 


W’aCfa) 


_L5Pa 

2/xa  dz 


Oa  -  ?a) 


(67) 


\  PdesP,,  ( t  3'des) 

Wdes(Jdes)  —  Tdes  (  ^des  0  ) 

Pdes  v  7  7 

The  air  side  pressure  drop  and  the  solution  thickness  are  calculated  from 
aPa  _  3  /y»a 

_  2Tfpa?3 


(68) 


(69) 


hies 


3/»dcs  Udes 

Pdesd  _ 


(70) 


The  results  illustrated  that  the  dehumidification  and  cooling  process  increased  with  the 
following  parameters:  low  air  Reynolds  number,  high  Cu  particle  volume  fraction, 
increases  in  the  height  and  length  of  the  channel,  and  decreases  in  the  channel  width.  Also, 
it  was  found  that  an  increase  in  Cu  volume  fraction  led  to  more  stability  of  the  solution. 

Ali  and  Vafai  [22]  presented  the  mathematical  and  numerical  formulations  to  evaluate 
the  effects  of  the  inclination  angel  of  parallel-and-counter-flow  on  heat  and  mass  transfer 
between  air  and  falling  desiccant  film,  with  Cu  nanoparticles  suspended.  As  Ali  et  al.  [21], 
the  same  assumptions  were  held  and  two  approach  models,  introduced  by  Xuan  and 
Roetzel  [1 1]  where  used  in  this  literature.  Then,  the  simulated  conditions  used  in  this  study 
are  summarized  as  follow:  low  and  high  air  Reynolds  number  for  the  inclined  parallel  flow 
arrangement,  and  low  and  high  air  Reynolds  number  for  the  inclined  counter  flow 
configuration.  The  problem  configuration  for  this  study  was  shown  schematically  in  Fig.  8. 

Similar  to  Ali  et  al.  [21],  the  general  governing  equations  (continuity,  momentum, 
energy,  and  mass  diffusion  for  air  and  the  desiccant  film)  for  both  parallel  flow  and  counter 
flow  can  be  given  by 


—  For  the  air 
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Fig.  8.  The  geometric  configurations  in 
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the  study  of  Ali  and  Vafai  [22]. 
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In  order  to  solve  the  governing  equations  of  the  air,  the  coordinate  transformations  are 
needed  for  both  inclined  parallel  and  counter  flow.  The  simulation  conditions  used  in  this 
work  are  in  the  range  of  low  and  high  air  Reynolds  number. 

From  the  numerical  results  it  could  be  concluded  that  the  dehumidification  and 
cooling  process  of  air  and  regeneration  process  of  the  liquid  desiccant  are  enhanced  with 
increasing  inclination  angle  for  both  parallel-and-counter  flow-channels.  Moreover,  the 
effective  thermal  conductivity  of  the  liquid  desiccant  varied  with  the  Cu  nanoparticles  and 
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the  thermal  dispersion  and  low  desiccant  concentration  led  to  enhancement  of  the 
dehumidification  process. 

Xuan  et  al.  [23]  proposed  the  thermal  Lattice  Boltzmann  model  to  display  the  flow 
features  and  heat  transfer  process  of  Cu-water  nanofluid  flowing  inside  a  channel.  The 
important  advantage  of  this  method  was  that  it  took  the  molecular  dynamics  into  account 
and  bridged  the  gap  of  microscopic  or  macroscopic  phenomena  of  the  nanofluids.  For 
application  of  this  method  to  model  the  thermal  fluid  flow  problem,  there  are  two  existing 
methods:  the  multispeed  (MS)  and  the  double-distribution-function  (DDF),  the  latter 
method  being  used  in  this  article.  This  led  to  two  separate  distribution  functions  as  follow: 
the  temperature  is  treated  as  a  passive  diffusing  scalar,  and  it  is  simulated  by  a  distribution 
function  independent  of  density  distribution,  respectively.  The  feature  point  of  the  Lattice 
Boltzmann  method  is  the  assumption  that  the  particles  are  mesoscopically  located  at  a 
series  of  lattices  and  their  distributions  correspond  to  the  Boltzmann  distribution.  For 
simulation  of  heat  transfer  characteristics  of  the  nanofluid,  the  viscosity  dissipation  term  in 
the  Lattice  Boltzmann  equation  for  energy  transport  will  be  neglected.  The  equation  is 
given  by 

\Ti(x,t)-  T.q(x,t)\ 

Tj(x  +  e,At,  t  +  At)  -  T,(x,  t)  =  -  1  ' - 1  ;  /  =  1,2, 3,4  (80) 

XT 

where  ex  is  the  lattice  velocity  vector,  tt  is  the  collision-relaxation  time  constant  for  the 
temperature  distribution.  Teq  is  the  equilibrium  function  for  the  temperature  which  is 
defined  as 

If  =  AtT  +  i  =  1,2, 3, 4  (81) 

cz 

where  c  is  the  reference  lattice  velocity  (Ax/ At) 

By  using  the  Lattice  Boltzmann  method  to  model  the  nanofluid,  the  dynamic  factors 
influencing  the  flow  behavior  of  the  nanofluid  must  be  considered  first.  Without  external 
forces,  the  nanoparticles  are  under  influence  of  the  buoyancy  force,  gravitational  force, 
Brownian  force,  drag  force  and  dispersion  force  resulting  from  the  repulsive  potential  [13]. 
Hence,  the  sum  of  the  total  forces  vector  acting  on  the  nanoparticle  per  unit  lattice  volume 
( V)  is  given  by 

F  =  n(F  \  +  Fb+Fd  +  Fh)/ V  (82) 

where  n  is  the  number  of  the  nanoparticles  in  a  given  lattice,  h\  is  the  dispersion  force,  /’B 
is  the  Brownian  force,  FD  is  the  drag  force,  and  Fu  represents  the  buoyancy  and 
gravitational  force,  respectively.  These  several  variables  are  defined  as 
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Fd  —  —6n[iaAu 


(87) 


(88) 


where  a  is  the  radius  of  the  nanoparticle,  A p  and  Ah  are  the  density  and  the  velocity 
difference  between  the  nanoparticle  and  the  base  fluid,  respectively,  X  is  the  friction 
coefficient,  kB  is  the  Boltzmann  constant,  G',  is  the  zero  mean  (unit  variance  independent 
Gaussian  random  number),  L  is  the  distance  between  center-to-center  of  the  particles,  and 
A  represents  Hamaker  constant  which  accounts  for  the  material  properties  but  not  for  the 
configuration  shape. 

The  numerical  results  showed  that,  for  1%  volume  fraction  the  Nusselt  number  of 
Cu-water  nanofluid  was  27%  higher  than  that  of  pure  water  at  the  same  Reynolds 
number.  Moreover,  fluctuations  of  the  Nusselt  number  of  the  nanofluid  taking  place  along 
the  main  flow  direction  may  cause  random  motion  of  the  particles  suspended  in  the  fluid. 


4.  Conclusion 

Evidently,  nanofluids  have  great  potential  for  heat  transfer  enhancement  and  are  highly 
suited  to  application  in  practical  heat  transfer  processes.  This  offers  an  opportunity  for 
engineers  to  develop  highly  compact  and  effective  heat  transfer  equipment.  Several 
published  articles  show  that  the  heat  transfer  coefficient  of  nanofluids  are  much  higher 
than  the  common-base  fluid  and  give  little  or  no  penalty  in  pressure  drop.  The  main 
reasons  for  the  heat  transfer  enhancement  of  the  nanofluids  may  roughly  be  listed  as 
follow:  the  suspended  nanoparticles  increase  the  thermal  conductivity  of  the  fluids,  and  the 
chaotic  movement  of  ultrafine  particles  increases  fluctuation  and  turbulence  of  the  fluids 
which  accelerates  the  energy  exchange  process.  Furthermore,  that  the  numerous 
correlations  for  single-phase  fluid  have  clearly  failed  to  predict  the  heat  transfer 
coefficients  of  nanofluids  may  be  caused  by  the  complicated  phenomena  that  coexist  in 
the  main  flow.  Flowever,  some  sophisticated  correlations  for  predicting  the  heat  transfer 
performance  of  nanofluids  have  been  established.  At  the  present  time,  the  uses  of  the 
nanofluids  ar  still  remaining  in  early  stages  of  development.  Therefore,  urgent  theoretical 
and  experimental  research  works  is  needed  in  order  to  clearly  understand  and  accurately 
predict  their  hydrodynamic  and  thermal  characteristics.  Generally,  many  researchers 
indicated  that  nanofluids  behave  like  pure  fluids  because  the  suspended  particles  are 
ultrafine.  Flowever,  at  present,  no  formulated  advanced  theory  exists  to  explain  the 
behaviour  of  nanofluids  by  considering  them  is  as  multi-component  materials. 
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